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Abstract 

Recent astro-physical observations indicate that the proton-electron mass ratio and the fine struc- 
ture constant have gone through nontrivial time evolution. We discuss their time variation in the 
context of a dilaton runaway scenario with gauge coupling unification at the string scale M s . We 
show that the choice of adjustable parameters allows them to fit the same order magnitude of both 
variations and their (opposite) signs in such a scenario. 

PACS numbers: 98.80.Cq 
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I. INTRODUCTION 



In unified theories of fundamental interactions, a variety of fundamental constants are not 
necessarily "constant" but can vary as a function of spacetime. Therefore, many experiments 
and observations have been done to test the constancy of various fundamental constants [l[ • 

Among them, several groups report nonvanishing time variation of some of the funda- 
mental constants. For example, Murphy et al. report time variability of the fine structure 
constant a by use of absorption systems in the spectra of distant quasars j2|. They found 
that the fine structure constant a was smaller in the past, 

— = (-0.543 ±0.116) x 10" 5 , (1) 
a 

for the redshift range 0.2 < z < 3.7, 1 though similar observations of other groups do not 
necessarily reproduce this result [H, 0] • The linear interpolation of such a change yields the 
rate of the change, 2 

- = C(10~ 15 ) yr' 1 = C(10- 65 )M G , (2) 
a 

where the dot represents the time derivative and Mq = 1/ y/8irG is the reduced Planck scale. 

The observations of H 2 spectral lines in the Q 0347-383 and Q 0405-443 quasars also 
suggest a fractional change in the proton-electron mass ratio \i = m p /m c , 

— = (2.4 ±0.6) x 10" 5 , (3) 

for a weighted fit 0, which implies that it has decreased over the last 12 Gyr. The linear 
interpolation of such a change yields the rate of the change, 

- ~ -2.0 x 10~ 15 yr^ 1 = -1.7 x 10~ 65 M G . (4) 

Thus, though there are still large uncertainties, the hints of the time variation of fundamental 
constants are found. 

On the other hand, from the theoretical point of view, it is natural to allow time and 
space dependence of fundamental constants. In fact, superstring theory, which is expected 
to unify all fundamental interactions, predicts the existence of a scalar partner <fi (called 
dilaton) of the tensor graviton, whose expectation value determines the string coupling 
constant g s = e 9 ^ 2 [?|]. The couplings of the dilaton to matter induces the violation of 
the equivalence principle and hence generates deviations from general relativity. Therefore, 
though the dilaton is predicted to be massless at tree level, it is usually assumed that it 
acquires a sufficiently large mass, associated with supersymmetry breaking, to satisfy the 
present experimental constraints on the equivalence principle. 

However, Damour and Polyakov proposed another possibility which can naturally rec- 
oncile a massless dilaton with experimental constraints They pointed out that full 
string-loop effects modify the four dimensional effective low-energy action as 



s = |^v^(^i?±^[2a0-(V0) 2 ]-^H0)^ 2 -^±---) • (5) 

1 In Q , the data is updated but the time variation has almost the same value Aa/a = (— 0.57±0.11) x 10~ 5 . 
In this paper, we use the units of c = h = 1. 
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Here -Bj(0) (i = g,(f>, F, . . .) are dependent coupling functions. String scale M s is given by 
M s = a ~ x l 2 . In the weak coupling limit g s — > (0 — > — oo), they are expanded into 

5,(0) = e -« + c f + c « e * + c« e 2 ^ + ■ • • , (6) 

which comes from genus expansion of string theory with B t = S n (n = 0, 1, 2, . . .). 

Assuming a universality of the dilaton coupling functions, that is, the existence of a value 
4> m of which extremizes all the coupling functions it has been shown that, during 

a primordial inflationary stage, the dilaton evolves towards the special value m , at which 
it decouples from matter (so-called "Least Coupling Principle") jlj. Subsequent (slight) 
change of the dilaton induces the time variation of fundamental constants. Note that the 
dilaton becomes almost homogeneous in space during inflation so that spatial variations of 
fundamental constants are expected to be much smaller than their time variations. 

On the other hand, in the infinite bare coupling limit a* — > oo (0 — * +oo), it is suggested 
that all the coupling functions have smooth finite limits [9[], 

5,(0) = C t + Ae'* = Ci(l + fa-*), (7) 



with di = Aj/Cj. Since A { (di) is expected to be positive in the "large N"-type toy model 
of 0], we assume that all Aj's are positive in this paper. All the coupling functions are 
extremized (minimized) at m = +oo. In this case, also, the dilaton evolves towards its 



fixed point at infinity during inflation so that it decouples from matter |10l . Illl ]. In Ref. 



Ill ], assuming the dilaton coupling to dark matter and/or dark energy, the magnitude of 
the time variation of the fine structure constant is estimated. 

In this paper, we discuss time variations of the proton-electron mass ratio and the fine 
structure constant in the context of a dilaton runaway scenario with gauge coupling unifi- 
cation at the string scale M s . We show that our model can account for the putative time 
variation of these constants of the same magnitude with the opposite signs. 

The rest of the paper is organized as follows. In the next section, we calculate time varia- 
tion of the proton-electron mass ratio and that of the fine structure constant by taking into 
account effects associated with thresholds in renormalization group running and variations 
in the vacuum expectation value (VEV) of the Higgs field. 3 We then apply it to the specific 
case of a runaway dilaton scenario in §111. In §IV we show that our model can not only fit the 
observations well but also is testable by the experiments to verify the equivalence principle. 
Finally §V is devoted to the conclusion. 



II. RENORMALIZATION GROUP ANALYSIS OF THE TIME VARIATION OF 
THE FUNDAMENTAL CONSTANTS 

In this section, we calculate time variation rates of fundamental constants such as the 
proton-electron mass ratio and the fine structure constant from a fundamental point of view 
using renormalization group analysis. As for the particle contents, we concentrate on the 
standard model and its minimal supersymmetric extension, although our discussions could 
proceed in the same way for more extended models as well. 



For related works in the context of GUT, see 12 1 
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A. (i 



First of all, we focus on the time variation of the proton-electron mass ratio /x = m p /m c . 
The time variation of \i is given by 

ji = m p rh e 
[i m p m e 

Though the proton mass m p depends not only on the QCD scale Aqcd but also on the 
masses of the up quark and the down quark, we set m p to be proportional to A QC d because 
these quark masses are much smaller than Aq C d- Then assuming |m u |, |rhd| <C |Aqcd|, the 
time variation of m p is given by 4 

rhp _ Aqcd , q n 
m p Aqcd ' 

The QCD scale Aqcd can be extracted from the Landau pole of the renormalization 
group equations as 

1 1 bi . ( M a \ 6 3 . /M 



^lnU4^ +^ln' msuSY 



«3(Aqcd) oix(M a ) 2n \M S vsy J 27r \ m t 

2tt \m h J 2tt \mj 2tt VAqcd 7 1 ' 

Here the parameters b{ are given by bi = (bi, b 2 , b 3 ) = (41/10, —19/6, —7), b^ = (b\, b%, b? 3 ) = 
(33/5,1,-3), and bl~ h = -23/3, b*~ c = -25/3, b% = -9. a x = a x {M s ) is the gauge 
coupling unified at the string scale M s , Msusy is the supersymmetry (SUSY) breaking scale, 
and m t ,m^,m c are the masses of top, bottom, and charm quarks, respectively. Reduction 
to the case of nonsupersymmetric theory would be obvious, that is, we take Msusy = M s . 
Then, the time variation of the QCD scale is given by 

^3 Aqcd _ _«x _^|_Ms h - b s 3 M S usy h m t 



27tAq C d Oi\ 2nM s 2vr M SU sy 2nm t 

| b% m c b\~ h ( rh b m t \ b\~ c f m c m h \ 
2tx m c 2ir \mb m t J 2ix \m c J 

Here and hereafter, for simplicity, we assume the universality of the time dependence of 
fermion masses, that is, the time variation of all the fundamental fermion masses except 
three light quarks, u, d, and s, 5 is identical, which is denoted by mf/ntf. Later we present a 
set of sufficient conditions to realize such a universality in the context of a dilaton runaway 
scenario. Under such a universality, the last two terms of the right hand side of the above 
equation are dropped. Then, the time variation of the proton mass is given by 

m p _ Aqcd _ 2tt a x 1 M s 4 Msusy 2 m f 

m p ~A Q cD~9ax 3M S 9 M SUS y 9m f ' 1 j 
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Note that the light quark masses may contribute to the proton mass, which may lead to the slight change 



of the coefficient in Eq.(T2D [13 1. 



The current masses of these light quarks are irrelevant to our analysis. 
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After all, the time variation of the proton-electron mass ratio \x = m p /m e becomes 



(jl _ 2-k d x 1 M s 4 M SUSY 7 rhf 

li ~ 9 a\ 3 M s 9 Msusy 9 m f ' 1 j 

where we have used the universality of the time variation of fermion masses m e /m c = rhf/mf. 



B. a 

Next, we discuss the time variation of the fine structure constant a. The renormalization 
group equations yield the scale dependence of the gauge couplings a« {i = 1,2) 

1 ± + WJM + w^Y (14, 



ai(M EW ) a x 2vr \ m susy/ 2tt \ M EW 

where Mew — Mz is the electroweak symmetry breaking scale. Assuming SU (5) GUT, for 
example, the fine structure constant a can be related to the gauge couplings (z = 1,2) at 
any scale as 

1 5 1 

+ 



a (M EW ) 3ai(M EW ) «2(M EW ) 



3a x ix \ Msusy J 6ir V Mew J 

for Mew < Msusy- After the electroweak symmetry breaking, charged fields acquire their 
masses. Therefore, taking their mass thresholds into account, the fine structure constant a 
is given by 

i = _i_ + w^ +E & n (j^y (16) 

a a(M EW ) 2ir \ m t ) 2ir \m, lt J 

Here, the third term in the right hand side corresponds to fermion mass thresholds, where 



fi = t,b,---,e, b a = 32/3 and ba 1+1 denotes beta-function coefficients for fermion mass 
thresholds between fj and fj+i. 

Then, the time variation of the fine structure constant is given by 

a 80 Ax 50 M s 257 M gugY 7 M EW 116 m { 



a 



TJa\ 9ttM s 54tt Msusy 2vr M ew 27vrm f ' 



where we have used the universality of the time variation of fermion masses except three 
light quarks u, d, s, and m u ~ m& ~ m s ~ Aqcd because light quarks are confined and 
effectively have a mass of the order of Aqcd- 

We have assumed that the gauge couplings are unified to ax- The gauge coupling unifi- 
cation is consistent with experimental values on the gauge couplings within the framework 
of the minimal supersymmetric standard model, but not in the (non-SUSY) standard model. 
For the latter case, we need some corrections to gauge couplings at M s . Such corrections can 
appear from gauge kinetic functions, which depend on moduli fields other than the dilaton. 
We assume that moduli-dependent corrections to the gauge couplings do not vary while only 
the dilaton varies in the time range relevant to our analysis. 
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C. Universality of Time Variation of Fermion Masses 

The four dimensional effective low-energy action related to the generation of fermion 
masses is given by 

S D J d i x^[B H ^){D,H^D^H u + B H ^){D,H d )W^H A 

• /Vi, /:! //; :: ( v ;./7;//., • //^i i: 0:Ji//,i • //^ , (is) 

where Jp = 7 M -D M , represents the covariant derivative, and j is the generation index. 
Canonically normalizing all fields yield effective Yukawa couplings and fermion masses, 



rriv. 



B„. 



Vd- = Vd- — TUd 



VdtVd 



^JBq.B^Bh^ y/2 ' 

~ B V°i Vc^d , , 

Uplike and downlike fermions acquire masses m u . = y u .v u /\^2 and = Ud^d/ V2 through 
the Higgs mechanism, respectively. Here y\ is a Yukawa coupling constant and v u and v d 
are the VEV of the up-type and down-type Higgs fields, respectively. Note that v u (v d ) is 
replaced by the VEVs of the standard Higgs field v in the case of the nonsupersymmetric 
standard model with a single Higgs doublet. 

Then, the time variation of fermion masses is given by 

"is. = s + ^ 

m Ui y Ui v u ' 

"hi. = + ^ 

m d] y di v d ' 

The Ve- Vd , s 

— — = ^ + — . (20) 
m e . y e . } v d 

Thus, the universality of the time variation of fermion masses is realized, for example, if the 
following conditions are satisfied, 

= 2/dk = 2/ei = m Vu = Vd = V 

2/uj Vd k 2/ei Hi v u v d v' 

where j, k, 1 are generation indices and the last equality comes from v 2 = v 2 + v 2 . Note that 
the second condition is unnecessary in the case of the nonsupersymmetric minimal standard 
model. Here, we have implicitly assumed that the scale dependence of the Yukawa couplings 
is negligible, that is, the time dependence of the Yukawa couplings is dominated by their 
dilaton dependence. 
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Here, we comment on radiative corrections on Yukawa couplings. Renormalization group 
effects on Yukawa couplings are the same among quarks except the top quark. Thus, ra- 
diative corrections do not violate the universal time variation of Yukawa couplings among 
quarks except the top quark. Furthermore, the mass ratios of other quarks to the top quark 
do not change drastically between M s and Mz, i.e. rri{(M B ) / m t (M B ) ~ rrif (Mz)/m t (Mz). 
Hence, even including radiative corrections, the universal time variation of quark Yukawa 
couplings would be a reasonable assumption, and such corrections on the time variation on 
Aqcd would be sufficiently small. The same discussion holds true for the universal time 
variation of Yukawa couplings only among leptons. However, radiative corrections on quark 
Yukawa couplings are different from those on lepton Yukawa couplings, because of correc- 
tions from a^. Such difference would be estimated as Y q /Y q —Yi/Yi = aax, where \a\ < 0(1), 
and it has some effect on a /a Eq. (jXTJ) , but it can be neglected compared with the first term 
in Eq. (fTTj). 



D. Electroweak scale, Higgs VEV, and SUSY scale 



Before investigating the time dependence of the VEV of the Higgs fields v/v , we discuss 
the time dependence of the electroweak symmetry breaking scale Mew characterized by the 
gauge boson mass Mz, 




M z = -vV + S' 2 = vJit [- ai + a 2 ), (22) 



where «j = ^(Mew), and g and g' are S77(2)l and U{1)y gauge coupling constants respec- 
tively. Then, the time variation of the electroweak symmetry breaking scale Mew/Mew is 
given by 

M EW _ 3di + 5d 2 t> 
M EW ~~ 2(3«i + 5a 2 ) v' 
Inserting Eq. (jl4|) into this equation yields 



M EW ^ 3a? + ha 2 2 a x 99aJ + 25aj M s 45aJ + 125q| M susy v ^ v 
M EW ~ 2(3«i + 5a 2 ) «| ~ 207r(3ai + 5a 2 ) M s + 247r(3ai + 5a 2 ) M SUSY + v ~ v' ^ ' 

where we have used a.\ ~ 1/60 <C 1 and a 2 ~ 1/30 <C 1. 

We now investigate the time variation of the VEV of the Higgs fields, v/v. First, we 
consider the nonsupersymmetric minimal standard model. The Lagrangian density related 
to the standard Higgs field is expected to read 

= (D.HYD^H - A (h*h - j\ , 



(25) 



with H = ^BftH, A = B\\/B^, and v 2 = B^v 2 . Assuming that v is intrinsic and has no 
time dependence, the time variation of the VEV of the standard Higgs field is given by 

* = (26) 
v 2B n y ' 
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Next, we consider the SUSY model. The neutral components of up and down sector 
Higgs fields, h u and h d , have the following potential, 

V = m\\h d \ 2 + m\\h u \ 2 + b(h d h u + h.c.) + \{g 2 + g' 2 ) (\h d \ 2 - \h u \ 2 ) 2 , (27) 



with 



m i = m H d + Hh, m 2 = m H u + Vh, (28) 

where m 2 H are SUSY breaking scalar masses squared of Higgs fields, h Ujd , and fin is the 
super symmetric mass parameter. In addition, the parameter b is also SUSY breaking pa- 
rameter with mass dimension two, that is, the so-called &-term. 
By using the stationary conditions, dV/dh U)d = 0, we obtain 

\{9 2 + 9> 2 = -ml -m\ + ^|^K " ml), (29) 
-6(tan/3 + cot/?) =m\ + m\, (30) 

where tan/5 = v n /v&. Furthermore, for a moderate and/or large value of tan/3, i.e. tan 2 /5 > 
0(1), these equations reduce to 

\{g 2 + g' 2 )v 2 = -ml (31) 

-6tan/3 = m\ + m\. (32) 

That implies that if the time variation of mass parameters is the same, i.e. 

-T- In u 2 tt = ^- In m 2 TT = ^-\nm%r =— In 6, (33) 
dt ^ H dt Hu dt Hd dt K ' 

tan/3 does not vary, that is, v n /v n = Vd/v<i- The above assumption ( 1331) might be plausible 
for mass parameters at tree level, but m 2 H has a significant radiative correction due to stop 
mass, 

* m Hu ^m(M s /m f ), (34) 

where is SUSY breaking stop mass. Thus, in general, the value of tan j3 varies in time, 
and the time variations of v u /v n and Vd/^d are different. To take into account this aspect, we 
have to consider the situation that the time dependence of up-type quark masses are different 
from those of down-type quark masses and lepton masses, and we have to introduce another 
parameter to represent such difference. Such extension of our analysis is straightforward 
and would enlarge a favorable parameter space. (Note that because of Msusy/A^susy the 
SUSY model has more degrees of freedom than the non-SUSY standard model.) Similarly, 
the time variation of v also depends on those of several values, m 2 Hu , yt, m~ as well 
as the gauge couplings. To simplify our analysis, we use the same parameterization as the 
non-SUSY model Eq.flH. 

Now let us discuss the time dependence of the SUSY breaking scale Msusy- Although it 
strongly depends on the SUSY breaking model, we give one example based on the gaugino 
condensation and gravity mediation model. We consider a hidden sector, in which a hidden 
gauge group with a coupling blows up and hence gauginos A a condensate at some scale 
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M c , which breaks the SUSY. Then, repeating the same argument as the case of the QCD 
scale, the condensation scale M c is given through the RG flow by 



M c = M s e «hCM s) 6 h ? ( X a \ a ) = M c 3 0, (35) 

where &h is the beta-function coefficient which, for example, is given by = 3N C for the 
gauge group SU(N C ). If this breaking is transmitted to the visible sector through the 
gravitational interaction, the SUSY breaking scale Msusy is given by 



In this case, the time variation of the SUSY breaking scale is given by 



Msusy = 8ttGM c 3 = -^§-. (36) 



Msusy _ 3 M C _ 67r d h (M s ) JA G 

M S usy M c M g b h ag(M 8 ) M G 

III. TIME VARIATION IN A RUNAWAY DILATON SCENARIO 
A. Runaway Dilaton 

Now, we estimate the time variation of the proton-electron mass ratio and the fine struc- 
ture constant in the context of a runaway dilaton scenario. From the four dimensional 
effective low-energy action ([5]), we have the following relations, 

B rW = (38) 

Since the coupling functions are given in Eq. ([7]), the dilaton dependence of the gravitational 
coupling G and the unified gauge coupling ax is given by 

M G = (87TG)- 1 = 2M*B g ((j)) = 2M S 2 C 9 (1 + d g e^), 

a" 1 = 8nB F ((p) = 8nC F (l + dpe^), (39) 

which leads to 

Mq 2(1 + ^) M s ~ 2 gC ^ + M S ' 

^ = 87rA F e~^0 > 0, (40) 
«x 

where we set > without generality and we have assumed <C 1. We regard the string 
scale M s as fundamental and hence it has no time dependence M s = 0. 

In the context of the dilaton runaway scenario, the sufficient conditions for the universality 
of the time variation of Yukawa couplings, (1211) . are satisfied, for example, in the case that 
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the dilaton dependent functions have the following properties, 

£? Qj = B Lk = B D = C D + = C D (1 + d D e^), 

5 Uj = 5 dk = B Cl = B S = C S + Ase-* = C s (l + 4e" ), 

B yu . = B yik = B yci = B y = C y + A y e~* = C y (l + d y e^), 

£? Hu = B Ud = B n = C H + A H e-^ = C H (1 + d&r*). (41) 

Hereafter we assume that the dilaton dependent functions satisfy the above conditions (|4TI) . 
In this case, the universal time variation of fermion masses rhf/rrif is given by 

fhf m v 

— = — + -, 42 

m f y f v 

where the universal time variation of Yukawa couplings yt is estimated as 

^~-tZe-*0, (43) 

with 

d = dy - (d D + rf s + d u )/2, (44) 

and ( 1261) reads 

^ = - 2(l + ^) ^-2 dHe V (45) 

Note that <i can take either a positive value or a negative one, which is an important point 
to account for the decline of \i. 

After all, the universal time variation of fermion masses is given by 

TL-t + i^-fn.*)^. (46) 

m f y { v V 2 / 
B. /i and d in the Runaway Dilaton Scenario 

Finally, we show that the choice of adjustable parameters allows them to fit the observed 
time variation of the proton-electron mass ratio and the fine structure constant in the dilaton 
runaway scenario, and give constraints on the parameters. 

In the case of the non-SUSY model (Msusy = 0), the time variation of the proton-electron 
mass ratio /x = m p /m e reads 

(i 2n ax 7 rh{ 
H 9 a\ 9 mf 

2vr . 



— e 
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(47) 



where we have used Eqs. ( 1401 and (1461) in the second equality As is seen in Eq. (144|) . d can 
be negative. On the other hand, the time variation of the fine structure constant is given by 

a 80 ax 7 M EW 116 rhf 
«2 ~~ 27 aJ ~ 27 M EW + 277 m? 
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80 d x 
27 of 
80 
27 £ 



116 rhf 




(48) 



As mentioned in the introduction, observations indicate a/a > and fi/fi < 0, which impose 
constraints on the parameters, 



29 



1607T 2 



, 43 , 

a H aw 

464 H 



< A 



P < 



7 



167T 2 



(49) 



These constraints can be easily satisfied if d < 0, |d| ^> <i H and < 29|<i|/(1607r 2 ). Though 
the time variations of the proton-electron mass ratio and the fine structure constant depend 
on the evolution of the dilaton field 0, from Eqs. ([I]), ©, fT4"7]) . and (14"8"|) . we can expect that 
observed variation can be fitted by taking our model parameters appropriately around order 
of unity. 

In the case of the SUSY model, the time variation of the proton-electron mass ratio 
fi = m p /m c becomes 



2tt a x 4 Msusy 



9 a\ 
9 



9 Msusy 
8ttA f ( 



7rhf 
9 nif 

- 12 



-4,, 



A F b h 



7 

2^ 



d 



du_ 
2 



Ad, 
7 



(50) 



where we have used Eqs. ( )37l) . (1401) and (1461) . On the other hand, time variation of the fine 
structure constant is given by 



a 



a- 



80 d x 257 Msusy 
27 a\ 54vr Msusy 

^e'U 8itA f ( 1 



7 Mew H6 rhf 

277T fflf 

29 



27tM E w 
771 A h 



80 A F 6 h y 20vr 



257 

d + ^-Ah - ^d g 



43 
464' 



(51) 



In this case, too, from Eqs. (pQ), ([3]), (l5"Uj) . and (1511) . we can expect that our parametrization 
can fit the observed variation naturally. In fact, the SUSY model has more degrees of 
freedom than the non-SUSY standard model because of Msusy/Msusy- Furthermore, if we 
introduce the time variation of the ratio of v u /v u to Vd/vd as well, we have a wider parameter 
space to fit the observed variation. 



IV. COMPARISON WITH OBSERVATIONS 

Having formulated the time variation of the fundamental constants in the particle-physics 
context and given their explicit form in the dilaton runaway scenario, we now solve cosmolog- 
ical evolution of the dilaton field to show our model can account for the observed variation. 
Before giving an explicit result, however, we must consider other experimental consequences 
of the dilaton coupling which impose a stringent constraints on the parameter space. 
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A. Experimental Constraints on Dilaton Coupling 



The dilaton coupling to hadronic matter induces deviations from general relativity: post- 
Newtonian deviations from general relativity and the violations of the equivalence principle 
[sl EH • After integration by parts, the action of the dilaton is rewritten as 



,5' 



d xy g 



BM ~ z{<p) 



a' 



d g e 



(V0) 2 - V(</>) 



(52) 



where Z{<j>) =C+- A^ and B g {<j>) = C g (l 

From the effective mass m{4>) of a test particle (composed of hadronic matter) in the 
Einstein frame metric = B g ((f))g^, using E q.(T2"]) . the strength of the coupling of the 
dilaton to hadronic matter, aw, is given by BE] 



C*had 



l 2B g (<p) d\nm(<p) 

z\<j>) d~4> 



' IBM 



din m p ((j)) 



ld\nB g (<j>) 
2 dd> 



16tt' 

~1T 



\d n 



(53) 



where the non-SUSY standard model is assumed in the second line. Since the mass of a test 
particle depends on the dilaton, the test particle experiences an acceleration, — Vlnm, in 
addition to the usual free-fall acceleration g, which results in the violation of the universality 
of free-fall, oihad is related to the Eddington parameter 7 measuring a post-Newtonian 
deviation from general relativity and to the Eotvos ratio r\ measuring the difference in 
accelerations, a/, between the two test masses (I = A,B) BE!: 
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1 



a 



'^ a hady 



had 



a A — a B 

a A + a B 



5.2 x 10 



-5^2 



a 



had 



-2.6 x 1(T 5 (7 - 1). 



(54) 
(55) 



The present experimental constraints are 7 — 1 = (2.1 ± 2.3) x 10 [14j and, rj = (—1.9 ± 



2.5) x 10~ 12 for A =Be and B =Cu [l5j, and 77 = (0.1 ±4.4) x 10~ 13 from the measurement of 
the composite-dependent differential acceleration of Earth and Moon toward the Sun, which 
are made of several materials such as Fe, Ni, Si, Met, Ollf 



B. Dilaton Evolution 



We numerically solve the cosmological evolution of the dilaton to give a concrete example 
which can account for the time variations of /1 and a. Its dynamics is determined by 
specifying the dilaton potential V(<fi). We assume the following form of V(<p) in accord with 
Eq. ([7]): V{4>) = Vo(l + d v e~^). We further assume that the magnitude of V is similar to 
that of the cosmological constant, Vq ~ 10~ 120 Mq, which is miniscule in the unit of the 
string scale in order for the dilaton to play the role of dark energy. 6 



6 We expect that the value of Vb may be determined by the expectation value of other fields \ which 
are orthogonal to the dilaton and the smallness of Vb may be related to the largeness of such fields as, 
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FIG. 1: The time evolution of \x and a (solid curve) in runaway dilaton scenario. The upper panel 
is fj, determined from the spectral lines of hydrogen molecules. The data are taken from [g|. The 
lower panel is a determined from quasar absorption lines at several redshifts. Open circles are the 
binned data from 0] (hi each bin the value of Aa/a is the weighted mean), crosses are from ji| 
and filled circles are from 0. The datum at z ~ 0.16,0.44 is the Oklo bound and the 187 Re 



bound 



20J, respectively. 



Before proceeding to numerical calculations, let us make analytic estimate. From ([53 
and (I55p . we find 

-I 2 



rj ~ I0~ 5 a 2 had ~ 10"V 



-2p 



2C, 



9 9\ 2 



1 , 

-Clr, 



<io 



-13 



(56) 



Assuming the last two factors are of order of unity, we find e 2< ^ < 10 7 which implies 0^8. 
Therefore the dilaton must have run away much beyond the Planck scale by now. In order 



Vo ~ M s 4 e _x . Alternatively, one may argue the smallness of Vb is due to a nontrivial vacuum structure 
which may lead to Vb = M^e~ SE where Se is a Euclidean action connecting two degenerate vacua 17 1 
Since the purpose of this paper is to demonstrate that it is possible to fit the observational results, we 
will not pursue this issue (the problem of the cosmological constant) any further here. 
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to reproduce the observed time variation, must change appropriately in the cosmological 
time scale. Specifically, from <^ and (1481 or (|5"T1) . we require a/a ~ ae~^<p ~ 10~ 65 , in the 
Planck unit with Mq = 1- Using the slow-roll equation of motion we estimate 



3# x/K 



V^e -0 ~ lO- bO 4e- , (57) 



in the same unit. From these results we find that d v must take a fairly large value, d v > 10 4 . 
We also find a similar value of d v from (j4j) . 

The presence of nonzero hardly affects the evolution of the dilaton since it is moving 
slowly: it only slightly facilitates the rolling of because the equation of motion of is 
divided by — A^e~^. However, the presence of nonzero d g greatly affects the evolution 
since it induces a negative effective potential for 0: Vr{4>) = —Cgdge'^R/a'. Vr practically 
vanishes during the radiation dominated epoch when R ~ —(a'/2B g )T ~ with T being 
the trace of the energy momentum tensor, and hence the dilaton does not move. It becomes 
non-negligible during matter dominated epoch. The mag nitude of Vr can be larger than 



V((j)) and hence can decrease rather than increase |18j], quite the opposite to what we 
want. Therefore, d g must be small enough in order to be consistent with the observations. 

Provided these conditions are satisfied, there are a wide allowed region in the parameter 
space that can account for the observed variation, so that we only give one specific example 
in Fig. [H where the time variations of fi (upper graph) and a (lower graph) are given. The 
data of /j, are taken from 0. The open circles in a data which are binned and averaged by 
redshift are taken from [2], crosses from 0] and filled circles from j^]. We also plot the bound 
of Aa/a from the analysis of isotope abundances in the Oklo natural reactor operated 2Gyr 
ago (corresponding to z ~ 0.16). We adopt a conservative bound by Damour and Dyson 
[191 ] . The bound on Aa/a over the age of the solar system ~ 4.6Gyr [z ~ 0.44) from the 
comparison of the 187 Re meteoritic measurements of the time averaged 187 Re decay rate with 
the laboratory measurements is also included. We adopt a conservative bound by Fujii and 
Iwamoto (Sj. 7 

In this example, we have considered a non-SUSY theory and taken d g = 0.01, = A^ = 
1, d v = 10 4 for dilaton couplings to gravitational part and Ap = 0.001, d = —0.4, dn = 0.01 
for those to gauge/matter part. C g is used to normalize the gravitational constant and set to 
be 0.5, which corresponds to M s ~ MqJ \/1C g ~ 2.4 x 10 18 GeV in this example. However, 
one should notice that the action (1521) is invariant under the following scale transformation: 
M s — > aM s ,(j) — > 0, C g — > Cg/c?,C^ — > C^/a 2 ,A^ A^/a 2 . Therefore, the dynamics of 
the dilaton is unchanged under the transformation. If we take a ~ 1/10, the string scale 
can be lowered to the GUT scale. The initial condition of the dilaton is <fi = 8, = at 
z = 10 10 . As is seen in the figure, this example can fit the observational data. Damour 
and Dyson obtained Aa/a = (—0.9 ~ 1.2) x 10~ 7 at z ~ 0.16 from the Oklo mine samples 
"ig| | . In our example, Aa/a = —7.8 x 10 -8 at z = 0.16, which satisfies the Oklo bound. 
The 187 Re bound is |Aa/a| < 1.5 x 10~ 6 at 2 ~ 0.44 0], but it is to be noted that the 



7 In these nulear bounds on varying a, if quark masses are allowed to vary, there may be correlations 
between the effects of varying a and those of varying quark masses, since nuclear effects depends both on 
electromagnetic forces and on nuclear forces. These effects would lead to different bounds on a. However, 
we will not consider such a possibility here for simplicity. 
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bound depends on the way of time dependence of the 187 Re decay rate [20|. 8 The bound 
is also satisfied in our example: Aa/a = —1.8 x 10~ 7 at z = 0.44. The comparison of 
different atomic clocks over several years yiel ds a bound on the present-day time variability 
of a (a/a = (—0.3 ± 2.0) x 10~ 15 yr _1 [23j), which is also satisfied in our example. We 
also note that this example reproduces the cosmic expansion law of the present universe 
correctly since the dilaton potential plays the role of the dark energy with = 0.74 and 
the equation-of-state parameter w ~ —0.91. 

Although the allowed region in the parameter space is large we mention that there is 
a tension to account for both the smallness of rj and the observed values of a/ a and \ij \i 
simultaneously, for the latter requires a relatively large <fi sourced by a large d v . In the 
present example, we find 7 — 1 ~ —2.9 x 10~ 9 and 77 ~ 7.6 x 10 -14 . 

Finally we comment on the time variation of the gravitational constant G. It is given by 
G/G = loihad&hadl (1 + a 1ad) < 0- However, its magnitude is very small due to the smallness 
of ahad- \GIG\ ~ — 2a/ M( i|d;/24a 2 | ~ 10 _18 yr~ 1 and safely satisfies the present experimental 
constraint [24J: G/G — (4 ±9) x 10 _13 yr _1 . The tests of the weak equivalence principle put 
severe limits on dilaton models. More precise experiments of the weak equivalence principle 
could discover its violation by a runaway dilaton. 



V. CONCLUSION 

Motivated by the observational evidence that indicates nontrivial time evolution of the 
fundamental constants such as the proton-electron mass ratio and the fine structure constant, 
we have theoretically calculated their time variations based on the standard particle-physics 
model and its supersymmetric extension with the help of the renormalization group ap- 
proach. We have employed several simple assumptions such as the universality of the time 
variation of the masses of the fundamental fermions as a first step. We have applied our 
formalism to a specific scenario of the runaway dilaton and found that we can account for 
the observed time variation with some natural choice of model parameters. 

We have also found that there is some tension between the observed magnitude of the 
time variation and the experimental constraint imposed by the verification of the equivalence 
principle. Indeed we typically find 77 to be larger than 10 -13 unless we adopt a sufficiently 
small value of Ap. One can regard this feature of our model as a prediction, that is, if one 
performs a more precise experiment on the equivalence principle to measure 77 with a higher 
accuracy, one would be able to discover its violation by a runaway dilaton. 
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